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This  report  is  «  sunsury  of  the  research  conducted  by  tbe 
principal  investigator.  Jay  I.  Valton,  under  A70SK  Grant  77-3JSP0 
for  the  years  4/1/77  to  3/31/82.  Much  of  the  work  described  in 
this  report  was  in  collaboration  with  Dr.  Arja  Machmaa  who  wee  a 
co-principal  Investigator  on  the  grant  froe  4/1/77  until  3/31/79. 
A  variety  of  theoretical  questions  concerning  contact  and  fracture 
phenoaena  for  elastic  and  viscoelastic  aaterial  and  associated 
mathematical  techniques  were  addressed  under  the  grant.  The 
particular  topics  studied  were: 

1.  quasi-static  contact  problens  for  power-law  linearly 
viscoelastic  material; 

2.  the  quasi-static  propagation  of  a  crack  In  a  power-law 
linearly  viscoelastic  material; 


3. 

4. 

5. 

6. 

A  brief  description  of  these  projects  and  a  summery  of  the 
research  accosqpllshaents  achieved  on  each  fellows  is  Che  went 
More  detailed  discussions  of  the  above  problem 


systems  of  generalised  Abel  integral  equations,  and 
simultaneous  dual  integral  equations; 

the  application  of  conservation  integral  techniques 
to  static  elastic  lnterfadel  edge  cracks; 

thermo-elastic  interfacial  cracks; 

•  ..  / 

the  dynamic,  steady-state  fracture  of  general  Unfairly  p*  •'/ 
viscoelastic  material.  '*<7 


■  ■  "  ' ,  ‘ 

various  annual  sad  technical  reports  that  have  been 


this  grant.  The  final  section  of  this  report  is  a 
of  abstracts  from  the  several  journal  publications  that  jeeps' 
the  research  supported  by  this  grant. 


SKTXQN  2.  saggy  of  awwyste  tetfrisr 


Initial  research  effort 


this  AfOSI  grant  invol 


study  of  the  quad-static  alldng  oontact  friction  bat**#*  two 


oolif  bodies,  one  of  which  is 


to  be. 


Hie  other  of 


which  is  aasuaad  to  be  nodellad  as  a  linearly  dacoalaatlc  solid. 


This  class  of  problaso  baa  significant  ongiaaarlng  appllaattoa  .fee 
the  prediction  of  the  frictional  force  between  two  bodies  in 
sliding  contact  where  one  exhibits  aechanical  hysteresis  (visco¬ 
elasticity)  and  the  other  is  very  ouch  harder  and  non  rigid.  Deo 
examples  of  such  Interactions  are  tires  (rubber)  against  paveaent 
(rock)  and  polymeric  materials  in  contact  with  rigid  objects  (as 
occurs  in  many  structural  applications  of,  say,  gr aphl te-epoxy 
composite  materials). 

Sliding  contact  problems  in  linear  viscoelasticity  theory 
have  received  considerable  attention  by  analysts  during  tbs  Last 
two  decades.  The  primary  objective  of  these  studies  ha* 
predict  the  dependence  of  the  contact  pressure  distribution, 
friction  coefficient  (or  mschsnlrsl  dissipation)  and 


on  tha  applied  feres,  the 


properties  of  the  other 


of  the  slider  and 


t  analytical  solution  found. 


which  the  neterlel  it  Harms  sited  ee  e 


i.e.,  the  vioeoeleetle 


nsotlel  function. of  tine  end  feieesn*e  t 


of  model  is  for  toe  si^le  to 


It  has  been  observed  experimentally  that  aaajr  viacoai— tCfcv 


materials  can  be  characCarlsad  ratbar  wall  with  a  a lap la  power* 
lav  aodal,  l.a.,  by  assuming  a  constant  Poisson's  ratio  and  a 
sbaar  aodnlua  that  la  a  power  function  of  tl«a.  Consequently ,  It 
was  daaaad  Important  to  attempt  an  analyala  of  tha  problaa  of  a 
rigid  lndaator  alldlng  ovar  a  power- law  linearly  viscoelastic 
half "plana .  It  vaa  alao  considered  to  ba  of  both  mathematical  and 
physical  intaraat  to  lacorporata  into  tha  nodal  a  Coaled  fric¬ 
tional  forca  between  tha  lndaator  and  tha  half -plana,  l.a. ,  under  the 
alidar  tha  ahearlng  tractlona  ara  proportional  to  tha  noxual  tractions 

A  clo«4d  form  solution  was  obtained  to  this  problaa  by  reduction 
to  a  ganaralised  Abal  type  Integral  aquation  for  tha  unknown  neraal 
traction  under  tha  slider.  This  equation  was  than  solved  by  a 
transformation  to  a  Rieeana-Hllbert  boundary  value  problen.  All 
Important  physical  parameters  ware  calculated  for  three  canonical 
lndaator  shapes:  a  parabolic  punch  (as  an  approximation  to  a 
cylinder),  a  wedge  punch  and  a  flat  punch.  This  work  appears  In 
tha  paper:  'The  sliding  of  a  rigid  lndentor  ovar  a  power  Inc 
viscoelastic  half-space , ”  J.  k.  Walton,  A.  If  schema  sad  ft*  A.  tBhapury, 
Quart.  J.  Mach.  Appl.  Math.,  (31)  #3,  197*. 

Subsequently,  the  problaa  of  a  rigid  lodtmter  »TldtM*dji^  .  . 


viscoelastic  layer  of  finite 


Math— ticelly,  tha  layer  problem  la  math  meow  difficult* 

"  '  4-- i* 

ie  the  paper:  _  "Rm  jUMpg  gf  •  : : " ' 


sad  J .  *. 


J.  dmpl*  Mach.  (100)  #1,  117*. 


2.  Subsequent  to  the  work  on  friction,  attention  we*  U 
on  a  study  of  tho  fracture  of  viscoelastic  aatarlal .  Tax 


aiailar  to  thoaa  sap: loyed  la  tt»  gassi-static  contact  analysis 
wax*  uoad  to  obtain  an  analytical  eolation  to  tho  problem  of  '^JSSsHraBr 


tha  propagation  by  a  rigid  wadga  of  a  2-dlaaaaioaal  crack  in  an 
inflaita  power-law  linearly  vlacoalaatic  notarial.  Ala  work 
appears  in  the  paper:  "Tha  propagation  of  a  crack  by  a  rigid  m 
in  an  infinite  power-law  viscoelastic  body,"  J.  1.  Walton  and 
A.  Nachman,  J.  Appl.  Mach.,  (46)  #3,  1979. 

The  content  and  significance  of  this  study  ia  perhaps  most 
succinctly  summarized  in  the  following  excerpt  from  the  cited 
paper. 

This  particular  problem  is  an  example  of  a 
physically  realizable  mode  of  fracture  to  which  the 
correspondence  principle  does  not  apply.  Moreover, 
the  close-dorm  solutions  obtained  here  provide  a 
method  for  experimentally  determining  physical  para¬ 
meters,  such  as  the  stress-intensity  factor, 
important  for  an  understanding  of  fracture  phenomena 
This  problem  has  been  studied  extensively  for 
elastic  material,  a  detailed  description  of  which 
appears  in  the  article  by  Barenblatt  (2).  Tbs  wedg¬ 
ing  problem  was  considered  by  Barenblatt  for  two 
reasons.  In  general,  Barenblatt  addressed  hlmealf 
to  the  Issues  of  unrealistic  elaguler  crack  tip 
stresses  and  Indeterminant  crack  lengths  smsswif#f#a 
in  fracture  in  classical  linear  elasticity,  ftp 
wedge  problem  exhibits  both  features  evsm  v^ ajh 
erack  is  "semi-infinite"  since  tbs  it** 

•bead  of  tha  wedga  la  finita.  la  paetiea 
thruat  of  Barenblatt* a  article  ia  that  1*  i 
hood  of  tha  creek  tip,  the  separated  *wrf*tit  'w!:  • ’  % 
experience  cohesive  nolecuier  forces 
a  cusp  profile  sad  a  finita  straps  field. 
accounts  for  tkis  by  As  pheitneaanlfiglnit 
of  assuming  that  for  equilibrium  crsdfcs  »  , 

normal  traction  acta  over  a  smell  interval  tufrfal  a*/  *  II 
■  crack  tip  with  a  constant  stress-intensity  fas***  -  *v  -  $$ 
chsrsctariatic  of  tbs  asterisl.  As  aat*d  :b^'l|i^^|y,-$i 
blatt,  tbs  finite  crack  length  is 


"■  j .  V, .  • j  .j  c  » u  /,  ■  1 


sag? •. iii 

„  t ,  r  '  W I  • ' 

L-'  ■'  • . 

i  xxmxm* 


is  deroralmad  SnSitfr 

"*  1  ■  ■  ^  ^  ^x,i 


,  ’  f'f V'v^lS 


vy*.1;  .'t-  - 


knowing  this  itrui-lnttuit;  factor  and  la  lade-  ,  ‘  . 

pendent  of  the  actual  distribution  of  cohesive 
forces.  The  utility  of  tbleobeervatlon  lias  In 
reversing  ths  Implication,  learn  ly,  If  sxpsrlnsn tally  a 
wedge  la  driven  into  an  clastic  notarial  and  ths 
lengtk  of  ths  free  surface  of  ths  resultant  equilibrium 
crack  is  measured,  then  ths  stress-intensity  factor 
■ay  be  easily  calculated. 

Xn  this  paper,  us  provided  simple  closed-fora 
expressions  for  ths  stress  and  diaplacarant  fields 
for  the  quasi-static  wedging  of  s  power  las  visco¬ 
elastic  material  under  the  Barenblstt  hypothesis.  This 
analysis  culnlnatss  in  s  rsletloa  among  the  notorial 
parameters,  the  material  stress-intensity  factor  and 
the  length  of  the  free  crack  surface. 

3.  Concurrent  with  the  viscoelastic  fracture  effort,  e  study  of 

a  broad  class  of  simultaneous  systems  of  generalised  Abel  integral 

equations  was  Initiated  and  completed.  This  study  was  motivated 

in  part  by  the  appearance  of  generalised  Abel  equations  in  the  power- 

law  viscoelastic  contact  and  fracture  work  and  in  part  by  the  fact 


that  such  systems  were  known  to  arise  in  the  analysis  of  certain 
simultaneous  systems  of  dual  integral  equations.  These  systems  of 
dual  Integral  equations  occupy  s  central  role  in  the  analysis  of 
elastic  bl-msdla  (Interfacial)  cracks.  This  work  is  contained  in 


a  series  of  two  papers:  "Systems  of  generalised  Abel  equations , 


M.  Lowengrub  end  J.  K.  Unites,  SIAM  J.  Math.  Anal.  (10)#*,  1979$ 
"Systems  of  generalised  Abel  integral  agrarians  with  Ippl IrtBltSSi 
to  simultaneous  dual  rale tiara,"  J.  t.  Belton,  ibid. 


*.  further  questions  on 


Specifically,  the  problem  of  dsteralnieg  the  i 
(a  quantity  necessary  for  a  prediction  of  the 
of  a  cracked  solid)  for  elastic  intarf acini  c 


WSm 


The  aseoelated  boundary  rales  problem  is  racy  «lf  fienlt  ge 


explicitly. 


by  wee  of  the 


by  L.  B.  f round  la  the  calculation  of  stress-intensity 

for  cracks  in  homogeneous  clastic  aedla),  closed  fora  expreseiaf;.:'  <5! 

<  -  ■  s 

for  the  dasirod  energy  release  rate  were  obtained.  The  technique  ‘ 

was  adapted  to  both  a  aode  III  and  a  nixed  sodas  I  and  II  analysis .  '% 

Part  of  this  work  was  coapletad  while  the  principal  investigator, 

J.  R.  Walton,  was  a  Visiting  Associate  Professor  with  the  Maths- 

aatics  Kssearch  Canter  at  the  University  of  Wisconsin,  sad  appears 

in  the  MIC  Technical  Syanary  Report  #2017.  Moreover,  tip  project 

resulted  in  the  two  papers:  "Energy  release  rate  calculations  for 

Interface  edge  cracks  based  on  a  conservation  integral,  M  A.  Her  fail 

and  J.  R.  Walton,  Int.  J.  Solids.  Struct.  (16),  I960;  "Energy 

release  rate  calculations  for  an  interface  Mode  III  edge  crack  based 

on  a  conservation  Integral,"  A.  Madmen,  J.  Tweed  end  J.  R.  Walton, 

Lett.  Appl.  Englng.  Sci.,  (19),  1981. 


5.  In  collaboration  with  Professor  M.  Lowengrub  of  the  Department 
of  Mathematics,  Indiana  University,  a  study  of  the mo-elastic 


interfacial  cracks  was  initiated.  Specifically,  work  is  in  progress 
and  nearly  c Depleted  on  the  determination  of  the  combined  theme- 
mechanical  stresses  induced  by  a  penny-shaped  crack  at  tip  inWrfsob 
of  two  bonded  dlaslailax  elastic  half -spaces,  the  crank  IS  fo*«*d 
open  by  a  ranotaly  ispossd  fat  link.  Such  a  study  has  tfottant 
Will XiTl Til  iMlUnttw  to  m  aMlmlt  lit 

_  a  .  .  ..  ...  _  .  .  .  .  .. .  .  .  .  ..  _  .  ..-  ...  -■  .1  .  '  ~J,-|  ;  nil 


turns  involving  conpoeita  namriala  wndsi 
Ibis  week  will  appear  in  a  paper  ft.  la 
sLittti  for  publications. 


..  V.  ‘  ’,i  TMk  i 
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imsi 


vV"-' 


f.  Ae  part  of. «  collaborative  research  effort  with 

;'^!gpp 

Schapery  of  tho  Pspsrtmsate  ofAsrospece  mi  Civil  - 

Texas  AAM  University,  a  study  of  tho  dynamic  fracture  of  gansrsl 
linearly  viscoelastic  matsrlal  woo  initiator.  Konovtr,  dynamic 

■•:  A'  -* 

viocoolootic  frocturo  is  now  tho  principal  focus  of  research  undor 


tho  current  AFOSR  grant  of  the  principal  investigator,  J.  R.  Helton 
Three  separate  dynamic  crack  problems  have  been  considered;  each 


la  described  below 


A.  An  analysis  of  the  anti-plane  strain  problea  of  a  steadily 


propagating,  seal-infinite  crack  in  an  Infinite  linearly  visco¬ 


elastic  body  which  is  subjected  to  a  distribution  of  shearing 


tractions  moving  with  the  crack  was  completed.  This  work  is 


contained  in  the  paper:  “On  the  steady  state  propagation  of  an 


anti-plane  shear  crack  in  an  infinite  general  linearly  viscoelastic 


J.  k.  Walton,  Quart.  Appl.  Math.,  April,  19S2.  The  aignl 


flcance  of  these  results  is 


Scientific  leports. 

B.  An  analysis  was  completed  of  the  dyneaic  Mode  Ill  crack  problem 
for  a  layer  of  finite  thickness  analogous  to  that  deserted  if- V! 

paragraph  A  for  an  infinite  body.  A  dosed  form  empryasisn  Sse 

*  .  O/;.- 

the  screes  intensity  factor  was"  obtained  for  very  gsmam|l  wiS'ilS^  ':, 
elastic  shear  modulus..  The  solution  and  analysis  are  dmch  mime .  ■  j 

complicated  for  e  layer  than  for  an  infinite  body,  hit  the 
■Albert  method  that  mss  enployod  yields,  s  form  for 
intensity  factor  that  illaminatee  the  eomhimed  ef fecto  ef  msCdrtld 
inertia  ami  boundary  intaractione-  A  report  cm.  this 
being  written  for  p#lleetlim«  , ,<"• 


,  * 
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C.  A  study  of  the  steady-state,  dynamic  Mode  I  crack  problem  to 
an  Infinite  general  linearly  vlecoelaatlc  body  analogous  to  the 
Mods  III  problsm  discussed  in  paragraph  A  was  Initiated.  The 
Mods  I  model  Is  more  Important  to  engineering  applications  than 
the  Mode  III,  but  introduces  significant  additional  mathematical 
complications.  Only  preliminary  results  have  been  obtained  to  date 
and  these  are  described  in  the  renewal  proposal  of  the  principal 
investigator  that  was  submitted  in  December,  1981. 
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SECT1 OH  3.  Lift  of  Abstract*. 


Walton,  J.  R. ,  A.  Nachman  and  R.  A.  Schapery,  "The  sliding  of  a 
rigid  indentor  over  a  power- lew  viscoelastic  half-space," 
Quart.  J.  Mach.  Appl.  Math.,  (31)  #3,  1978. 

Suasury — Closed  fora  solutions  are  obtained  for  the  problem 
of  a  rigid  asperity  sliding  with  Coulomb  friction  over  a 
power- law  viscoelastic  half-space.  The  dual  integral  equa¬ 
tions  relating  the  unknown  normal  traction  under  the  contact 
Interval  (also  unknown)  to  the  unknown  normal  displacement 
outside  the  contact  interval  are  solved  by  first  reducing 
the  system  to  a  generalized  Abel  integral  equation  and  then 
appealing  to  the  theory  of  Riemann-Hilbert  boundary-value 
problems.  The  physical  quantities  of  interest  (eg.  the  co¬ 
efficient  of  sliding  friction)  are  determined  for  the  three 
canonical  indentors:  a  parabolic  punch,  a  wedge  punch,  and 
a  flat  punch.  The  analysis  predicts  singularities  in  the 
normal  stress  field  for  certain  power-law  materials  even  for 
the  smooth  parabolic  indentor. 

Nachman,  A.  and  J.  R.  Walton,  "The  sliding  of  a  rigid  indentor 

over  power-law  viscoelastic  layer,"  J.  Appl.  Mech.  (100)  #1, 
1978. 

Summary — The  problem  of  the  sliding  of  a  rigid  asperity  over 
a  power  law  viscoelastic  layer  is  examined  in  the  realistic 
limit  of  infinite  (dimensionless )  layer  thickness,  For  a 
contact  interval  of  unit  length,  asymptotic  expansions  for 
the  normal  traction  over  the  interval  together  with  several 
other  physically  relevant  quantities  (e.g.,  the  friction 
coefficient)  are  developed  in  terms  of  an  appropriate  asymp¬ 
totic  sequence  of  powers  of  the  (dimensionless)  layer 
thickness. 

Walton,  J.  R.  and  A.  Nachman,  "The  propagation  of  a  crack .by  a 
rigid  wedge  in  an  infinite  power  law  viscoelastic  body,"  J. 
Appl.  Mech.,  (46)  #3,  1979. 

Summary — Closed-form  solutions  are  obtained  for  the  problem 
of  a  crack  propagated  by  a  2-dimensional  rigid  wedge  of 
finite  thickness  in  an  infinite  power  law  viscoelastic  body. 
The  distance  from  the  crack  tip  to  the  point  of  contact  with 
the  wedge  is  determined  from  die  assumption  that  in  the 
vicinity  of  the  crack  tip,  forces  of  cohesion  act  to  produce 
a  smoothly  closing  crack  with  continuous  stresses.  Finally, 
a  simple  formula  expressing  the  stress-intensity  factor  as  a 
function  of  the  speed  of  penetration  and  the  length  of  the 
free  crack  surfaee  is  exhibited. 


Lowengrub,  M.  and  J.  R.  Walton,  "Systems  of  generalized  Abel 
Equations , "  SIAM  J.  Math.  Anal.  (10)  #4,  1979. 

Abstract — Certain  mixed  boundary  value  problems  arising  in 
the  classical  theory  of  elasticity  lead  to  the  solution  of 
certain  systems  of  generalized  Abel  integral  equations.  A 
method  is  presented  where  these  systems  are  reduced  to 
uncoupled  pairs  of  Riemann  boundary  value  problems.  Closed 
form  solutions  are  obtained.  We  also  demonstrate  how 
general  systems  of  dual  relations  (given  in  terms  of  Erdelyi- 
Sneddon  operators  of  fractional  integration)  may  be  reduced 
to  these  systems  of  Abel  equations. 

Walton,  J .  R. ,  "Systems  of  generalized  Abel  integral  equations 
with  applications  to  simultaneous  dual  relations,”  SIAM  J. 
Math.  Anal.  (10)  #4,  1979. 

Abstract — A  method  is  presented  for  solving  certain  systems 
of  generalized  Abel  integral  equations  by  constructing  an 
equivalent  system  of  singular  integral  equations.  An  appli¬ 
cation  is  then  given  to  a  class  of  simultaneous  dual  relations 
of  a  type  arising  in  bimedla  fracture  problems  in  elasticity. 
The  equations  discussed  in  this  paper  generalize  those 
considered  in  an  earlier  paper  of  Lowengrub  and  Walton  [SIAM 
J.  Math.  Anal.,  this  issue,  pp.  794-807]. 

Nachman,  A.  and  J.  R.  Walton,  "Energy  release  rate  calculations 
for  Interface  edge  cracks  based  on  a  conservation  integral," 
Int.  J.  Solids.  Struct.  (16),  1980. 

Abstract — A  certain  conservation  integral,  the  so  called  M- 
lntegral,  recently  exploited  by  L.  B.  Freund  in  the  calcula¬ 
tion  of  stress  intensity  factors  for  cracks  in  homogeneous 
elastic  media,  is  applied  to  the  calculation  of  energy 
release  rates  for  interface  edge  cracks.  Specifically,  for 
an  edge  crack  along  the  interface  between  two  elastic  wedges 
of  different  opening  angles  and  dissimilar  elastic  properties, 
and  that  is  subjected  to  point  loads  at  the  apex,  a  relation 
is  derived  among  the  length  of  the  crack,  the  energy  release 
rate  of  the  crack,  the  applied  loads,  the  wedge  angles  and 
the  material  parameters. 

Nachman,  A.,  J.  Tweed  and  J.  Walton,  "Energy  release  rate  calcula¬ 
tions  for  an  lntsrfacs  Mods  III  edge  crack  based  on  a 
conservation  Integral,"  Lett.  Appl.  Engng.  Sci.,  (19).  1961. 

Abstract — The  M-lntegral  is  applied  to  the  calculation  of 
energy  release  rates  for  Interface  edge  creeks  of  the  Mode 
III  type.  Specifically,  for  an  edge  crack  along  the  inter** 
face  between  two  elastic  wedges  of  different  opening  angles 
and  dissimilar  elastic  properties,  and  that  is  subjected  to 
point  loads  at  the  apex,  a  relation  is  derived  among  the 
length  of  the  crack,  the  energy  release  rate  of  the  crack, 
the  applied  loads,  the  wedge  angles  and  the  material  para¬ 
meters. 


12 


Walton,  J.  R.,  "On  the  steady-state  propagation  of  an  anti-plane 
shear  crack  in  an  infinite  general  linearly  viscoelastic 
body,"  Quart.  Appl.  Hath.,  April,  1982. 

Abstract — The  steady-state  propagation  of  a  semi-infinite 
anti-plane  shear  crack  is  considered  for  a  general  infinite 
homogeneous  and  isotropic  linearly  viscoelastic  body. 

Inertial  terms  are  retained  and  the  only  restrictions  placed 
on  the  shear  modulus  are  that  it  be  positive,  continuous, 
decreasing  and  convex.  For  a  given  integrable  distribution 
of  shearing  tractions  travelling  with  the  crack,  a  simple 
closed-form  solution  is  obtained  for  the  stress  intensity 
factor  and  for  the  entire  stress  field  ahead  of  and  in  the 
plane  of  the  advancing  crack.  As  was  observed  previously 
for  the  standard  linear  solid,  the  separate  considerations 
of  two  distinct  cases,  defined  by  parameters  c  and  c*, 
arises  naturally  in  the  analysis.  Specifically,  c  and 
c*  denote  the  elastic  shear  wave  speeds  corresponding  to 
zero  and  infinite  time,  and  the  two  cases  are  (1)  0  <  v  <  c* 

and  (2)  c*  <  v  <  c,  where  v  is  the  speed  of  propagation 
of  the  crack.  For  case  (1)  it  is  shown  that  the  stress 
field  is  the  same  as  in  the  corresponding  elastic  problem 
and  is  hence  independent  of  v  and  all  material  properties, 
whereas,  for  case  (2)  the  stress  field  depends  on  both  v 
and  material  properties.  This  dependence  is  shown  to  be  of 
a  very  elementary  form  even  for  a  general  viscoelastic  shear 
modulus . 
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